Atomic four-wave mixing: fermions versus bosons 
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We compare the efficiency of four-wave mixing in quantum degenerate gases of bosonic and 
fermionic atoms. It is shown that matter-wave gratings formed from either bosonic or fermionic 
atoms can in principle exhibit nearly identical Bragg-scattering, i.e. four-wave mixing, properties. 
This implies that effects such as coherent matter-wave amplification and superradiance can occur 
in degenerate fermi gases. While in the boson case the Bragg resonance is clearly due to 'Bose 
enhancement', in the case of fermions the resonance is due to constructive many-particle quantum 
interference. 

PACS numbers: 03.75.Fi, 03.75-B 



The experimental realization of nonlinear atom op- 
tics [0 is one of many recent advances made possible 
by the achievement of atomic Bose-Einstein condensa- 
tion. The demonstration of atomic four- wave mixing , 
the discovery of BEG superradiance , and the devel- 
opment of coherent matter-wave amplification are 
all examples of nonlinear wave mixing involving atomic 
Bose-Einstein condensates. One proposed interpreta- 
tion of these phenomena, in which the generation of a 
new atomic and/or optical side- modes is attributed 
to 'Bose-enhancement', predicts that these effects would 
not be seen in Fermi systems. In contrast, a second in- 
terpretation which attributes the effect to Bragg scat- 
tering from atomic matter-wave gratings, predicts that 
four-wave mixing with fermions should, in fact, be pos- 
sible. 

In this letter we address this fundamental issue and 
compare the Bragg-scattering properties of quantum- 
degenerate bosonic and fermionic matter-wave gratings. 
Three issues need to be addressed: First, there is the 
question of whether it is indeed possible to create a high- 
quality matter-wave grating in a quantum degenerate 
Fermi gas. Second, one needs to determine the scat- 
tering properties from such a grating. Lastly, a precise 
quantum-mechanical interpretation of atomic four-wave 
mixing must be formulated. Our analysis shows that 
the establishment of a fermionic grating is indeed possi- 
ble, and furthermore leads to scattering properties prac- 
tically indistinguishable from those of a BEG with the 
same mean density profile. We then present a physi- 
cal interpretation of this result, contrasting the Bose en- 
hancement responsible for four-wave mixing in the boson 
case to a many-particle quantum interference effect in the 
fermion case. 

We consider a system consisting of two scalar fields 
with annihilation operators denoted as >f'i(r) and •^2(1"). 
The first field contains N identical bosonic or fermionic 
atoms from which a matter-wave grating will be formed. 
The second field contains a single test particle which will 
probe the scattering properties of the grating. This test 



particle might be an atom, in which case the scattering 
properties are related to atomic four-wave mixing exper- 
iments, or it could be a photon, in which case the results 
would be applicable to phenomena such as BEG superra- 
diance and matter-wave amplification. The commutation 
relations for the two fields are therefore given by 



i^i(r)#|(r') - ^i,,w]{Y'%{v) = S,,6\r - r') 
!^,(r).^,(r')-M„'^,(r')'^»(r) = 0, 



(1) 



where the coefficients fiij vary depending on the quantum 
statistics of the two fields. If both fields are bosonic we 
have fj,ij = 1. If, however, the grating is fermionic then 
either the test particle is an atom of the same species in 
a different internal state, in which case fj,ij = —1, or it is 
a different type of bosonic (fermionic) particle, in which 
case /ill = -1, M12 = /^2i = 1, and ^^22 = 

The two fields are subject to the free Hamiltonians Hi 
and 7^2, respectively, and are coupled via a two-body 
interaction of the form 



V = A / rfVif|(r)ifJ(r)ip'2(r)iZ'i(r), 



(2) 



which describes equally well atom-atom collisions in the 
s-wave scattering approximation or the effective inter- 
action between ground-state atoms and far off-resonant 
photons. For simplicity we neglect collisions between 
atoms in the grating. 

At present, the primary technique for establishing den- 
sity gratings in BEG is by use of a 'beam splitter' based 
on two-photon Bragg transitions between center-of-mass 
states . Each atom is thereby transferred into a coher- 
ent superposition of its initial state and a copy of that 
state displaced in momentum space by the two-photon 
recoil kick. Quantum interference between these two mo- 
mentum groups then results in a typical 'standing wave' 
density modulation. 

For the purpose of this letter we specialize to the case 
of gratings formed from harmonically confined gases at 
zero temperature. We first introduce the 3-dimensional 
harmonic oscillator states 
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/ (,x)4'l3{m){y / ^y)(l>',{m){z / 



(3) 



where is the oscihator length along the j-axis, a{m), 
(3{m), and 7(771) are the quantum numbers of the rrith 3-d 
oscillator energy level, and is the normalized rith 1-d 
harmonic oscillator energy level. We further define the 
creation operators al^Ck) for the momentum side-modes 
of the TTith oscillator state as 



«L(k) = J 



(4) 



At r = the states of A^-atom Bose and Fermi gases are 
[at (0)]^|0)/VAn and Ul=l ^1(0)10), respectively. Thus 



lm=0 ' 

immediately after the Bragg pulse is applied, the state of 
the matter-wave grating is given for bosons by 



N 



|V'B) = [2^iV!]-^/^[4(0) + e^'«S(K)J |0) 
and for fermions by 



(5) 



= [2^]-V2 -Q [st^(o) + e^<^al{K)] |0), (6) 

where 9 and hK. are the relative phase and momentum 
transfer imparted on the atoms by the Bragg splitter. 
Henceforth we assume 6 = for simplicity. We note that 
these states are only normalized in the limit that KK. is 
much larger than the rms spread in momentum of the ini- 
tial Bose and Fermi clouds, a condition which we assume 
to hold. It is this condition, in fact, which is sufficient to 
guarantee a high-quality grating at wavelength 2t:/K. 

To compare the Bragg-induced gratings in degenerate 
Bose and Fermi gases we compare mean atomic densities 
p(r) = (#|(r)#i(r)) for the two states (||) and (1). By 
making use of the commutation relations (|l]) we find 



/9B(r) =iV|^o(r)|'[l + cos(K-r) 



and 



PFir)=Y. l^m(r)P[l + cos(K.r)]. 



(7) 



(8) 



which shows that in both cases, the effect of the Bragg 
splitter is to superpose a density modulation [1 + cos(K • 
r)] to the mean density of the initial atomic cloud. Fig. 
1 shows the integrated density p{x, z) = J dy p{r) for 
the case of spherically symmetric harmonic traps with 
oscillator lengths for the Fermi gas and = ^/2Axf 
for the Bose gas, where Axf is the rms spread of the 
initial Fermi cloud. This choice of £,b guarantees that the 
rms widths of the two clouds are equal in r-space. The 
parameters for Fig.l are K = (20/^)2 and N = 86, and 
X and z are given in units of ^. Figure la shows the 
result for the BEC, while Fig. lb shows the result for 
the degenerate Fermi gas. 



This example establishes that it is possible to construct 
very similar matter-wave gratings from degenerate Bose 
and Fermi gases. It should be noted that while their 
mean densities are nearly identical in coordinate space, 
they are very different in momentum space. As the phase- 
space density of the BEC is much higher than that of the 
degenerate Fermi gas, the bosonic grating is significantly 
more localized in momentum space. 

The next step is to compare the scattering properties 
of these two gratings. Of particular interest here is the 
issue of Bose enhancement, which is oftentimes invoked as 
the cause of stimulated scattering in bosons. In order to 
address this question, we consider the case where a single 
test particle is incident on the matter-wave with some 
wave-vector ko and use perturbation theory to compute 
the probability P(k, t) that it is scattered into a given 
state k after some time t. The initial state of the system 
is therefore taken as 



where 



|V'(0))=ct(ko)|V'B.F), 



(9) 



(10) 



creates a plane-wave test particle of momentum kg in 
the quantization volume V . The probability at time t 
of finding the test particle in the plane-wave state k is 
therefore 

F(k,i) = (V(t)|ct(k)P^c(k)|V(0>, (11) 
where is the solution of the Schrodinger equation 

\ip{t)) = exp \-i{-Hi + n2 + v)t/n^ IV'(O)) (12) 

and Vn is the A^-atom projection operator 

Vn = f d^n...d^rN^l{rN)...^l{ri)\0) 



(0|#i(ri)...^i(rjv) 



(13) 



This scattering problem can be solved analytically by 
expanding the solution (12) to first order in the param- 
eter A of Eq. (H). This yields the scattering probabili- 
ties after the test particle as been scattered once by the 
grating. In order to proceed we first specify the free 
Hamiltonians Hi and Ti.2- Our approach is to assume 
that the states (^^(r) exp[jk • r] are approximate eigen- 
states of the first-quantized version of Tii , with energy 
i?,„(k) huim + h^k"^ /2M where fiujm is the energy 
of the TOth trap eigenstate and M is the atomic mass. 
These states are essentially low lying levels of the the 
harmonic trap shifted in momentum space by Kk.. For 
ultracold atoms it is reasonable to neglect the evolution 
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of the wavepackets for times small compared to the os- 
cillator length divided by the velocity hk/m. For these 
short times the overlap between the initial wavepacket 
and the state it evolves into remains approximately unity. 
In addition, we take the plane- waves (I/a/F) exp[ik • r] 
to be eigenstates of the first-quantized version of ii2 with 
energy Ti^k"^ /2m. 

An analytic expression for P(k, t) can be derived in a 
relatively straightforward manner. For k 7^ ko and |K — 
ko I K we find that the leading order (in A) contribution 
to Eq. (O) is given for the boson case by 



|Fi(k,i)|2 + |F2(k,t)|^ 



+ (iV-l)|i^i(k,i)+^^2(k,i)nGoo(k-ko)P 
+ (7V-l)|F2(k,i)|2|Goo(k-ko-K)|2] , (14) 

and for the fermion case by 



PF{k,t) « [|Fi(k,t)|2 + |F2(k,i)P 

Zn V 

+ \Fi{k,t)+F2{k,t)fg{k^ko) 
+ |F2(k,i)|2g(k-ko-K)] , 



where the time-dependent functions 

sin[2lr(k-ko)-kt] ^« 



Fi(k,t) 



[2lr(k-ko)-k] 



(k-ko)-kt 



(15) 



(16) 



and 



sin rTr-(k — kn) • (k — K)tl . r , z^\t 
F2(k, t) = r^^"' lJ.e«2;;r(k-ko)-(k-K)t 



[5|,(k-ko).(k-K)] 



(17) 



give the effects of energy conservation, and the functions 
G„„(k) = ^Jd^r ^;,(r)e^'^-"-)^„(r) (18) 



and 



^^(k) = — 

^ ' N 



^ ^ G,n,,i(k) 



N-1 



E \Gmnikf 



(19) 



describe the shapes of the Bragg resonances. The approx- 
imation indicates that we have dropped several (negligi- 
ble) terms from the exact expression. 

The first two terms in Eqs. ( [l^ ) and ( p^ ) correspond 
to spontaneous scattering, whereas the third and fourth 
terms correspond to 0th order (small angle) and 1st or- 
der Bragg resonances. In Figure 2 we plot the scattering 
probabilities in the k^-kz plane for the case kg = Kx. 
Figure 2a shows the relative probabilities for the bo- 
son matter-wave grating of Fig. la and Fig. 2b shows 



those of the fermion grating of Fig. lb. We have cho- 
sen t = .025(2m^2/?i) so that the Bragg resonance is 
slightly sharper than the energy conservation envelope. 
The large spot at /c^ = corresponds to small-angle (0th 
order) scattering, while the smaller spot at k^ — 20 to 
the 1st order Bragg resonance. The inset figures show 
close-ups of the 1st order resonances, i.e. the four-wave 
mixing signal. We see, perhaps surprisingly, that these 
two microscopically very different gratings produce al- 
most identical scattering cross-sections. This strongly 
suggests that (at least for large N) it is the mean atomic 
density alone which determines the efficiency of four- 
wave mixing processes. We note that for = 1 there is 
no Bragg-resonance, only spontaneous scattering. This 
shows that the Bragg resonance is a many particle effect, 
simply having a modulated mean density is not sufficient 
in the case of a single-atom grating. 

While both types of matter-wave gratings result in 
practically identical scattering distributions, the under- 
lying physical processes which produce the Bragg reso- 
nances seem at first glance to be quite different. We can 
gain considerable physical intuition about them by con- 
sidering the simple case TV = 2 in more detail. In this 
case, the states of the Bose and Fermi gratings are 



1 



2V2 



V2-^ [at(0)]'|0)+2at(0)aS(K)|0) 



V2^ [4(K)]'|0) 



(20) 



and 



1^/ 



1 r 



4(0)4(0)10) +aJ(0)aI(K)|0) 



4(K)aI(0)|0)+aS(K)aI(K)|0) 



(21) 



where we have written each term as an amplitude times 
a normalized state. In this representation the gratings 
can be thought of as quantum superpositions of different 
initial states. If the test particle is scattered into state 
k, then one of the atoms in the grating must absorb a 
recoil momentum kick of ?i(ko — k). This is accomplished 
by acting on the states ( pO| ) and ( pl| ) with the operator 
J d^r'Fl{r) exp[i(ko — k) • r]if'i(r). When the Bragg con- 
dition k = ko + K is satisfied, the resulting final states 
of the gratings are 
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at(-K)aS(0)|0)+aJ(-K)aS(K)|0) 



^/2i= [4(0)]'|0)+4(0)4(K)|0) 



(22) 



where the factor \/2 on the third term can be identified 
with 'Bose enhancement', and 



a;(-K)al(0)|0)-(-4(0)al(-K)|0) 



3 



4(-K)aI(K)|0)+aS(0)al(0)|0) 
4(0)at(0)|0)+aS(K)aI(-K)|0) 
4(0)aUK)|0)+at(K)at(0)|0) 



(23) 



In either case, the total transition probabihty is the sum 
of the squares of the amphtudes for all distinguishable 
"paths". If two final states are indistinguishable (e.g., 
the fourth and fifth terms of (p3|)) then the amplitudes 
must be added before taking the square modulus and 
adding to the other paths. 

In both cases, therefore, the peak relative probability is 
5/2, but in one case, it is attributed to Bose enhancement 
and in the other to constructive quantum interference. 
We observe that if the Bragg condition strongly violated, 
the Bose-enhancement factor disappears from IV'b)', and 
that term is replaced by two distinguishable states with 
weights equal to the other states. In the fermionic case, 
on the other hand, the two indistinguishable states be- 
come distinguishable. In both cases, however, the result- 
ing relative probability is reduced to 2. 

This example illustrates that the Bragg resonance in 
the case of a BEC can be attributed to Bose enhance- 
ment, whereas in the fermionic case it can be attributed 
to quantum interference between "paths" which lead to 
indistinguishable final states. That both mechanisms 
lead to precisely the same enhancement, a fact which 
is readily verified for = 3 and illustrated for = 86 
by Fig. 2, is due to the fact that Bose enhancement, 
when viewed in first-quantized many-particle quantum 
mechanics, is simply a constructive quantum interference 
where many initial states (the different terms under ex- 
change of particle labels) lead to a single final state (when 
all atoms are in same state there is only one term even 
after symmetrization) . Hence, effects which can be in- 
terpreted as Bragg-scattering from atomic-matter wave 
gratings, such as atomic four-wave mixing, BEC super- 
radiance, and matter-wave amplification, can work as ef- 
ficiently in principle in both Bose and Fermi systems. 
Because of the connection between Bose enhancement 
and quantum interference we find that, remarkably, it is 
possible to set up A'^-atom superposition states in Fermi 
systems which mimic the effects of Bose stimulation. 
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FIG. 1. The integrated mean density p(x,z) is shown for 
the case of degenerate Bose (Fig. la) and Fermi (Fig. lb) 
matter-wave gratings. The parameters are A'' = 86 and 
K = (20/5), S'lid the dimensionless variables x and ^ have 
been scaled to the oscillator length 5. 



FIG. 2. The relative scattering probabilities are shown for 
the case ko = Kx and ky = 0. Figures 2a and 2b correspond 
to boson and fermion matter-wave gratings, respectively. The 
inset figures show close-ups of the 1st order Bragg-resonances. 
The dimensionless parameters and kz have been scaled to 
1/e 
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This figure "figl.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0008264vl 
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